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Abstract 



Wc present a multiwavelength study of the NGC 281 complex which contains the young cluster IC 1590 
at the center, using deep wide-field optical U BVI C photometry, slitlcss spectroscopy along with archival 
data sets in the near-infrared (NIR) and X-ray. The extent of IC 1590 is estimated to be ~6.5 pc. The 
cluster region shows a relatively small amount of differential reddening. The majority of the identified 
young stellar objects (YSOs) are low mass PMS stars having age < 1 — 2 Myr and mass 0.5-3.5 M Q . The 
' slope (r) of the mass function for IC 1590, in the mass range 2 < M /Mq < 54, is found to be —1.11 ±0.15. 

Qh, The slope of the X-band luminosity function (0.37 ±0.07) is similar to the average value (~0.4) reported 

q for young clusters. The distribution of gas and dust obtained from the IRAS, CO and radio maps indicates 

clumpy structures around the central cluster. The radial distribution of the young stellar objects, their 
ages, A(H — K) NIR-excess, and the fraction of classical T Tauri stars suggest triggered star formation 
at the periphery of the cluster region. However, deeper optical, NIR and MIR observations are needed to 
have a conclusive view of star formation scenario in the region. The properties of the Class 0/1 and Class 
II sources detected by using the Spitzer mid-infrared observations indicate that a majority of the Class II 
sources are X-ray emitting stars, whereas X-ray emission is absent from the Class 0/1 sources. The spatial 
C ' distribution of Class 0/1 and Class II sources reveals the presence of three sub-clusters in the NGC 281 

West region. 

Key words: Galaxy: open clusters and associations: individual (IC 1590) - stars: formation - stars: 
luminosity function, mass function - stars:pre-main-sequcnce 



Introduction cess and makes an interesting and important topic of star 

formation. The formation of massive stars feeds energy 

H II regions have been studied quite extensively in re- back into the nearby environments, irradiating, heating 

cent years on account of their close association with star and compressing the remain of the natal molecular cloud. 

*_h . formation. There seems to be two modes of star forma- This feedback can have either destructive or constructive 

tion associated with H II regions depending on the initial effects, but it is not clear which dominates in a given cloud 

density distribution of the natal molecular cloud. One is or overall in a galaxy. Many examples exist in our Galaxy 

the cluster mode which gives birth to a rich open clus- as well as in other galaxies where spatial distributions of 

ters and the other is the dispersed mode which forms young stellar objects (YSOs) and their ages suggest for 

only loose clusters or aggregates of stars. Presumably, triggered star formation (see e.g. Walborn et al. 2002, 

the former takes place in centrally condensed, massive Oyc ct al. 2005, Deharveng et al. 2005, Sharma ct al. 

clouds, whereas the latter occurs in clumpy, dispersed 2007, Chauhan et al. 2009). 

clouds (see e.g., Ogura 2006). These clusters/aggregates The HII region NGC 281/Sh2-184 (a 20 oo = 00 ,l 52 m , 

of stars emerging from their natal clouds can be the lab- ^2000= ±56° 34' or 1=123°. 07, b= -6°. 31) is located at 

oratories to address some of the fundamental questions of a relatively high Galactic latitude and has the centrally 

star formation. Trends in their evolutionary states and located cluster IC 1590 (Guetter & Turner 1997, Henning 

spatial distribution can help distinguish between various et al. 1994). The brightest member of IC 1590 is an 

star formation scenarios such as spontaneous or triggered O-type Trapezium-like system HD 5005, whose compo- 

star formation. Triggered star formation is a complex pro- ncnt stars HD 5005ab (unresolved), HD 5005c, and HD 

5005d have spectral types of 06.5 V, 08 V, and 09 V, 
respectively (Walborn 1973; Abt 1986; Guetter & Turner 
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1997). Despite many measurements, the distance esti- 
mates of NGC 281 varies from 2.0 kpc to 3.7 kpc (cf. 
Sato et al. 2008). Recently, using the VLBI observations 
of an associated H2O maser source Sato et al. (2008) de- 
rived a trigonometric parallax of 0.355 ±0.030 milli arcsec, 
corresponding to a distance of 2.81 ±0.24 kpc. 

The NGC 281 region provides an excellent laboratory 
for studying in detail star formation through the inter- 
action of high mass stars with their surrounding cloud. 
Of special interest in this region is the possibility of trig- 
gered star formation occurring on two different scales: the 
large-scale (~300 pc) supernovae-triggcred formation of 
the first generation OB stars and their associated cluster 
(Megeath et al. 2002, 2003), and the subsequent, trig- 
gered sequential and ongoing star formation 011 a smaller 
scale (~1-10 pc); the latter is taking place in an adjoin- 
ing molecular cloud (NGC 281 West) probably through an 
interaction with an HII region (the NGC 281 nebula) ex- 
cited by the first generation OB stars (Elmegreen & Lada 
1978; Megeath & Wilson 1997; Sato et al. 2008). 

The southwestern quadrant of the NGC 281 nebula is 
obscured by the adjoining molecular cloud NGC 281 West. 
Ongoing star formation in NGC 281 West is indicated by 
the presence of H2O maser emission and IRAS sources 
within this cloud near its clumpy interface between the 
HII region. This star formation may have been triggered 
by the interaction of the molecular cloud with the HII re- 
gion (Elmegreen & Lada 1978; Megeath & Wilson 1997). 
The NGC 281 molecular cloud complex was mapped both 
in 12 CO (J=l-0) and 13 CO (J=l-0) emission lines by Lee 
& Jung (2003). The central radial velocity of the NGC 
281 West molecular cloud, Vlsr = 31 kms -1 (Lee & Jung 
2003) agrees well with that of the H2O maser emission in 
the cloud (Sato et al. 2007). Megeath et al. (2002, 2003) 
suggested that this cloud complex was formed in a frag- 
menting super-bubble, which gave birth to the first gener- 
ation OB stars, and these OB stars have then ionised the 
surrounding gas which subsequently triggered next gener- 
ation star formation in the neighboring clouds (Sato et al. 
2008). 

Though both low-mass and high-mass star-forming re- 
gions can be studied at variety of wavelengths ranging 
from radio waves to X-rays, however most of the present 
knowledge about the HII region/open cluster NGC 281 
has been inferred from studies outside the optical region. 
Henning et al. (1994) made a multiwavelength studies of 
the NGC 281/IC 1590 region including Stromgren pho- 
tometry of the bright cluster stars. The first detailed 
UBV CCD photometry of 279 stars for the cluster was 
published by Gucttcr & Turner (1997). Their photom- 
etry terminates at a magnitude limit that is marginally 
brighter than the expected brightness of pre-main se- 
quence (PMS) and T Tauri stars in the embedded clus- 
ter region. Keeping above discussion in mind we feel 
that NGC 281 is an appropriate target for a deep and 
wide field optical/infrared photometry. In this paper, we 
present deep wide-field optical UBVI C data and slitlcss 
spectroscopy. We supplement them with archival data col- 
lected from the surveys such as Chandra, Spitzer, 2MASS, 



IRAS and NVSS (NRAO VLA Sky Survey). Our aim is 
to understand the global scenario of star formation under 
the effects of massive stars in the whole NGC 281/IC 1590 
region. In Section 2, we describe our optical CCD photo- 
metric and slitless spectroscopic observations and briefly 
the data reduction. In Section 3, we discuss the archival 
data set used in the present study. In the ensuing sec- 
tions, we present the results and discuss star formation 
scenarios in the NGC 281 region. 

2. Observations and data reduction 

2. 1 . Optical photometry 

The CCD UBVI C observations of the NGC 281 region 
were obtained by using the 105-cm Schmidt telescope 
of the Kiso Observatory, Japan on November 21, 2004 
and November 27, 2005. The CCD camera used a SITe 
2048 x 2048 pixel 2 TK2048E chip having a pixel size of 
24 /im. At the Schmidt focus (f/3.1), each pixel corre- 
sponds to 1.5 arcsec and the entire chip covers a field of 
~ 50 x 50 arcmin 2 on the sky. The read-out noise and 
gain of the CCD are 23.2 e~ and 3.4 e~ /ADU, respec- 
tively. A number of short and deep exposure frames were 
taken. The average FWHM of star images was found to 
be ^3 arcsec. The observed region is shown in Fig. 1. The 
Kiso data were standardized by observing the cluster to- 
gether with the standard stars in the SA 98 field (Landolt 
1992) 011 07 January 2005 using the 2048 x 2048 pixel 2 
CCD camera mounted on the f/13 Cassegrain focus of the 
104-cm Sampurnanand telescope of Aryabhatta Research 
Institute of Observational Sciences (ARIES), Nainital. In 
this set up, each pixel of the CCD corresponds to ~ 0.37 
arcsec and the entire chip covers a field of ~ 13 x 13 
arcmin 2 on the sky. To improve the signal to noise ra- 
tio, this observation was carried out in the binning mode 
of 2 x 2 pixel. The read-out noise and gain of the CCD are 
5.3 er and 10 e~/ADU respectively. The FWHMs of the 
star images were ~ 2 arcsec. The log of these observations 
is given in Table 1. 

The CCD data frames were reduced by using computing 
facilities available at ARIES, Nainital. Initial processing 
of the data frames was done using the standard tasks avail- 
able from IRAF 1 and ESO-MIDAS 2 data reduction pack- 
ages. Photometry of cleaned frames was carried out by us- 
ing the DAOPHOT-II software (Stetson 1987). The PSF 
was obtained for each frame by using several uncontam- 
inated stars. Magnitudes obtained from different frames 
were averaged. When brighter stars were saturated on 
deep exposure frames, their magnitudes have been taken 
from short exposure frames. We used the DAOGROW 
program in constructing of an aperture growth curve re- 
quired for determining the difference between aperture 
and profile fitting magnitudes. Calibration of the instru- 
mental magnitudes to those in the standard system was 
done by using the procedures outlined by Stetson (1992). 

1 IRAF is distributed by National Optical Astronomy 
Observatories, USA 

2 ESO-MIDAS is developed and maintained by the European 
Southern Observatory. 
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To translate the instrumental magnitudes to the 
standard magnitudes the following calibration equations, 
derived using a least-square linear regression, were used: 

u = U + (7.004 ± 0.004) - (0.005 ± 0.006) (U - B) + (0.431 ± 
0.005)X, 

b = B + (4.742 ± 0.005) - (0.035 ± 0.004) (B - V) + (0.219 ± 
Q.004)X, 

v = V + (4.298 ± 0.002) - (0.038 ± 0.002) (V - I) + (0.128 ± 
0.002)X, 

i = I c + (4.701 ± 0.004) - (0.059 ± 0.003)(V - J) + (0.044 ± 
0.003)X 

where U,B,V and I c are the standard magnitudes and 
u,b,v and i are the instrumental aperture magnitudes nor- 
malized for 1 second of exposure time and X is the air- 
mass. We have ignored the second-order colour correc- 
tion terms as they are generally small in comparison to 
other errors present in the photometric data reduction. 
The standard deviations of the standardization residuals, 
A, between standard and transformed V magnitude and 
(U — B), (B — V) and iV — I) colours of the standard stars 
are 0.006, 0.025, 0.015 and 0.015 mag, respectively. Short 
exposure data of the cluster region taken on the stan- 
dardization nights, were standardized by using the above 
equations and coefficients. The standard magnitudes and 
colours of more than 50 stars obtained from these short 
exposures were further used to standardize the deep ob- 
servations taken with Kiso Schmidt. The standard devi- 
ations of the residual of secondary standards are of the 
order ~0.02 mag. The typical DAOPHOT errors in mag- 
nitude as a function of corresponding magnitude in dif- 
ferent pass-bands for the Kiso Schmidt observations are 
found to increase with the magnitude and become large 
(> 0.1 mag) for stars fainter than V ~ 20 mag. The mea- 
surements beyond this magnitude were not considered in 
the analysis. 

2.1.1. Completeness of the data 

To study luminosity functions (LFs)/ mass functions 
(MFs) , it is very important to make necessary corrections 
in data sample to take into account the incompleteness 
that may occur for various reasons (e.g. crowding of the 
stars). We used the ADDSTAR routine of DAOPHOT 
II to determine the completeness factor (CF). The pro- 
cedures have been outlined in detail in our earlier works 
(Pandey et al. 2001, 2005). Briefly, the method consists 
of randomly adding artificial stars of known magnitude 
and position into the original frame. The frames are re- 
reduced using the same procedure used for the original 
frame. The ratio of the number of stars recovered to those 
added in each magnitude interval gives the CF as a func- 
tion of magnitude. In the case of optical CCD photometry, 
the incompleteness of the data increases with magnitude 
as expected. The CF as a function of V magnitude is 
given in Table 2. Table 2 indicates that our optical data 
have a 95% completeness at V ~ 16.5 mag, which corre- 
sponds to a stellar mass of ~ 2 M Q for a PMS star having 
an age of ~2 Myr (cf. Fig. 11). 



2.1.2. Comparison with previous studies 

We have carried out a comparison of the present pho- 
tometric data with those available in the literature. The 
difference A (literature - present data) as a function of 
V magnitude is given in Table 3. The comparison indi- 
cates that the present V mag and B — V colour are in 
good agreement with the CCD and photoelectric photom- 
etry by Guetter & Turner (1997), whereas the A(U - B) 
shows a systematic variation with the V magnitude in the 
sense that the present (U — B) colours become blue with 
increasing V magnitude. 

2.2. Slitless grism spectroscopy 

Spectra of some PMS stars, specifically classical T- 
Tauri stars (CTTSs) show emission lines, among which 
usually Ha is the strongest. Therefore, Ha surveys have 
often been used to identify PMS stars. We observed the 
NGC 281 region in the slitless mode with a grism as 
the dispersing element using the Himalayan Faint Object 
Spectrograph Camera (HFOSC) instrument during two 
observing runs on 10 October 2005 and 16 August 2006. 
This yields panoramic images where the star images are 
replaced by their spectra. A combination of a 'wide Ha' 
interference filter (6300 - 6740 A) and Grism 5 (resolu- 
tion = 870) of HFOSC was used without any slit. The 
central 2K x 2K pixels of the 2K x 4K CCD were used 
in the observations. The pixel size is 15 /*m with an im- 
age scale of 0.297 arcsec pixel -1 . The observed sky area 
is shown in Fig. 1 as a large white box which was cov- 
ered by four field-of- views of ~ 10 x 10 arcmin 2 each. For 
each field-of- view we secured three spectroscopic frames of 
longer exposure with the grism in, and one direct frame of 
shorter exposure with the grism out for the identification 
purpose. The log of the observations is given in Table 1. 
Emission line stars with enhancement over the continuum 
at the Ha wavelength arc visually identified. 

3. Archival datasets 

3.1. 2MASS 

Near-infrared (NIR) JHK S data for point sources in 
the NGC 281 region have been obtained from the Two 
Micron All Sky Survey (2MASS) Point Source Catalogue. 
The 2MASS data reported to be 99% complete up to ~ 
16, 15, 14.7 mag in J,H,K$ bands respectively 3 . To secure 
the photometric accuracy, we used only the photometric 
data with the quality flag ph-qual=AAA, which endorses 
a S/N> 10 and photometric uncertainty < 0.10 mag. The 
NIR data are used to identify the Classical T-Tauri stars 
(CTTSs) and Weak line T-Tauri stars (WTTSs) (cf. §4.3). 

3.2. CHANDRA X-ray data 

Since YSOs are very strong X-ray emitters (as strong 
as log Lx/Ltoi ~ 10~ 3 ) and they can be detected behind 
column densities as large as Njii ~ 10 23 cm 2 (Linsky et al. 
2007), X-ray imaging of star forming regions and young 
clusters is valuable for identifying these sources. 

3 See http://www.ipac.caltech.edu/2mass/rclcases/allsky/doc/ 
sec6_5al.html 
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3.2.1. Observation 

Chandra observed the NGC 281 region on three occa- 
sions for 62.6 ks (Obs ID 5424, on 2005-11-10 @ 18:19:27 
UT), 23.5 ks (Obs ID 7206, 2005-11-08 @ 13:41:46) and 
13.1 ks (Obs ID 7205 @ 22:40:54 UT). The aim point of 
the array was a 200 o= 00 ?l 52 m 25 s .2, 5 2000 = +56°33'47".6, 
and the satellite roll angle (i.e. the orientation of the CCD 
array relative to the north-south direction) was 225°. 5 
for all observations. The exposures were obtained in the 
very faint data mode with a 3.2 s frame time using the 
ACIS-I imaging array as the primary detector. ACIS-I 
consists of four front illuminated 1024 x 1024 CCDs with 
the pixel size of ~ 0.492 arcsec and the combined field of 
view of « 16.9 x 16.9 arcmin 2 . The S2 and S3 CCDs in 
ACIS-S were also enabled, however in the present study 
we used the ACIS-I data only. The detailed information 
on Chandra and its instrumentation can be found in the 
Chandra Proposer's Guide (POG) 4 . To detect sources, we 
have merged the event-list files of all observations. The 
X-ray observed region is again shown in Fig. 1 by a black 
box. 

3.2.2. Data reduction and Source detection 

We analyzed the data reprocessed by the Chandra 
X-Ray Center on 2006 April 5 (ASCDSVER 7.6.7.1). 
The data were reduced by using the Chandra Interactive 
Analysis of Observations (CIAO; Fruscione et al. 2006) 
software (ver. 4.1; CALDB ver. 4.2). Light curves from 
the on-chip background regions were inspected for large 
background fluctuations that might have resulted from 
solar flares, and none were found. We have filtered the 
data for the energy band 0.5 to 7.5 keV. After filtering 
in energy, the time integrated background is 0.11 counts 
arcsec -2 . Source detection was performed on the merge- 
event list by using the 1.7 keV exposure map with the 
PWDetect 5 code (Damiani et al. 1997), a wavelet-based 
source detection algorithm. The significance threshold 
was set to 5 ct so as to ensure a maximum one spuri- 
ous source per field. We have detected 379 sources, out of 
which 9 sources either fell on the unexposed areas of the 
CCD or were doubly detected. This implies that a total 
of 370 X-ray sources were detected in the NGC 281 field. 
An IDL-based program ACIS Extract (AE; Broos et al 
2010) was used to extract the photons from each candi- 
date source in a polygonal region which closely matches 
with the local PSFs. The source free regions around the 
source were considered as background. AE provides the 
Poisson probability of not being a source. We have not 
considered those sources which have the probability of be- 
ing non-existence > 0.01. Sixteen such sources were found 
in the catalogue. After removing these sources the cata- 
logue consists of 354 X-ray sources. Further, the median- 
detcctcd photon energy for the point sources were- deter- 
mined by using the ACIS Extract software package. We 
estimated the background AGN rate within the Chandra 
field of view using the Chandra Deep Field (Brandt ct al. 
2001). At the 0.5-2 keV limiting flux of 5.6 x 10 -16 erg 

4 See http://asc.harvard.edu/proposcr/POG 

J See http://www.astropa.unipa.it/progctti_riccrca/PWDetect/ 



s -1 cm -2 , we expect to find 79 to 127 background objects 
within the Chandra field of view. The optical, 2MASS 
and IRAC counterparts of the X-ray sources were searched 
within a match radius of 1 arcsec and the data are given in 
Table 4. A sample of the table is given here, whereas the 
complete table is available in the electronic form only. Out 
of 354 X-ray sources, 193 and 90 sources have NIR and 
optical counterparts respectively. All the optical counter- 
parts of X-ray sources have NIR counterparts also. The 
location of X-ray sources in NIR colour-colour diagram 
has been used to identify the probable WTTSs/ Class 
III sources. The completeness of the X-ray data has not 
been estimated. Since we are using the WTTSs/ Class III 
sources to study the spatial distribution of these sources 
and to support the results obtained on the basis of rather 
complete optical and NIR data, the incompleteness of the 
X-ray data will not have any significant effect on the re- 
sults presented in this study. 

3.3. Spitzer IRAC data 

The Spitzer mid-infrared (MIR) surveys have enabled 
detailed censuses of YSOs in star forming regions. The 
classification of young stars as protostellar Class I or more 
evolved Class II sources with optically thick discs is best 
accomplished by using their broadband spectral energy 
distributions (SEDs) (Mucnch et al. 2007). 

We have used archived MIR data observed with Infrared 
Array Camera (IRAC) . We obtained basic calibrated data 
(BCD) using the software Leopard. The exposure time of 
each BCD was 10.4 sec and for each mosaic, 72 BCDs 
have been used. Mosaicking was performed by using 
the MOPEX software provided by Spitzer Science Center 
(SSC). All of our mosaics were built at the native instru- 
ment resolution of 1.2 arcsec pixel -1 with the standard 
BCDs. In order to avoid source confusion due to crowd- 
ing, PSF photometry for all the sources was carried out. 
We used the D AO PHOT package available with the IRAF 
photometry routine to detect sources and to perform pho- 
tometry in each IRAC band. The FWHM of every detec- 
tion is measured and all detections with a FWHM >3.6 
arcsec are considered resolved and removed. The detec- 
tions are also examined visually in each band to remove 
non-stellar objects and false detections. The sources with 
photometric uncertainties < 0.2 mag in each band were 
considered as good detections. A total of 347 sources were 
detected in the 3.6 and 4.5 \xm bands, whereas only 35 
sources could be detected in all the four bands. 

Aperture photometry for well isolated sources was done 
by using an aperture radius of 3.6 arcsec with a concen- 
tric sky annulus of the inner and outer radii of 3.6 and 
8.4 arcsec, respectively. We adopted the zero-point mag- 
nitudes for the standard aperture radius (12 arcsec) and 
background annulus of (12-22.4 arcsec) of 19.670, 18.921, 
16.855 and 17.394 in the 3.6, 4.5, 5.8 and 8.0 fxm bands, 
respectively. Aperture corrections were also made by us- 
ing the values described in IRAC Data Handbook (Reach 
et al. 2006). The necessary aperture correction for the 
PSF photometry was then calculated from the selected 
isolated sources and were applied to the PSF magnitudes 
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of all the sources. The 2MASS, optical and X-ray coun- 
terparts of the IRAC sources were searched for within a 
match radius of 1 arcsec. These counterparts are given 
in Table 5. A sample of the table is given here, whereas 
the entire table is available in the electronic form. The 
completeness of the data in the 3.6, 4.5, 5.8 and 8.0 fxm 
bands having S/N > 5 (error < 0.2 mag) is found to be ~ 
16.0, 15.5, 13.0 and 12.0 mag, respectively. 

34. IRAS 

The data from the IRAS survey in the four bands (12, 
25, 60 and 100 fun) for the NGC 281 region have been 
used to study the spatial distribution of warm and cold 
interstellar dust. One cold IRAS point sources is identified 
in the cluster region and its details are given in Table 6. 

4. Results 

4-1. Structure of the cluster 

4.. 1.1. Isodensity contours 

Internal interaction due to two-body encounters among 
member stars and external tidal forces due to the Galactic 
disc or giant molecular clouds can significantly influence 
the morphology of clusters. However, in the case of young 
clusters where dynamical relaxation is not important be- 
cause of their young age, the stellar distribution can be 
considered as the initial state of the cluster that should 
be governed by the star formation process in the parent 
molecular cloud (Chen et al. 2004). To study the mor- 
phology of the NGC 281 cluster, we plotted isodensity 
contours using the 2MASS data as well as table of the 
identified YSOs (cf. §4.3) in Fig. 2. The isodensity con- 
tours indicate an elongated morphology for the cluster. 
It is interesting to point out that sub-structures can be 
clearly seen in the 2MASS data (left panel of Fig. 2) 
towards the south-west of the cluster as well as in the 
south-east. 

4-1-2. Radial stellar surface density profile 

To find out the extent and radial stellar density profile 
of IC 1590 we used the 2MASS data of K < 14.3. First the 
cluster center was determined by using the stellar density 
distribution in a 100 pixel wide strip along both the X 
and Y directions around an initially eye estimated center. 
The point of the maximum density obtained by fitting 
the Gaussian distribution is considered as the center of 
cluster. It is found to be a (200 o) = 00 h 52 m 39 s .5 ± 1.0 s , 

5(2000) =+56°37 46"±15". 

To determine the radial surface density profile we as- 
sumed a spherical symmetry of stellar distribution and 
divided the cluster into a number of concentric circles. 
The projected radial stellar density in each concentric an- 
nulus was obtained by dividing the number of stars by 
its area. The densities thus obtained are plotted in Fig. 
3. The error bars are derived by assuming that the num- 
ber of stars in a concentric annulus follows the Poisson 
statistics. The horizontal dashed line in the plot indicates 
the density of contaminating field stars, which is obtained 
from the reference region ~ 15 arcmin away toward the 
northwest from the cluster center (01(2000) = 00 h 51 m 16 s .0, 
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5(2000) = +56°46 45 ). The extent of the cluster R c i is de- 
fined as the projected radius from the density peak to the 
point at which the radial density becomes constant and 
merges with the field star density. R c i from the optical 
data as well as from the NIR 2MASS data is estimated as 
~ 8 arcmin (^6.5 pc for a distance of 2.81 kpc). 

The observed radial density profile of the cluster was 
parametrized by following the approach of Kaluzny & 
Udalski (1992). The projected radial density profile p(r) 
is described as: 




where the core radius r c is the radial distance at which 
the value of p(r) becomes half of the central density po. 
The best fit obtained by the x 2 minimization technique 
is shown in Fig. 3. Within the uncertainties the model 
reproduces well the observed radial density profile of IC 
1590. The core radius r c comes out to be 1.7 ±0.4 arcmin 
(1.4 ±0.3 pc) and 2.0 ±0.3 arcmin (1.6 ±0.2 pc) for the 
optical and 2MASS data, respectively. 

4-2. Interstellar extinction 

4-2.1. Reddening 

The interstellar extinction in the cluster region is 
studied by using the (U — B)/(B — V) two-colour di- 
agram (TCD) shown in Fig. 4 where zero-age-main- 
sequence (ZAMS) from Schmidt-Kaler (1982) is shifted 
along the normal reddening vector having a slope of E(U — 
B)/E(B - V) = 0.72. The distribution of stars shows a 
small amount of differential reddening (E(B — V) ~ 0.2 
mag) in the region with the minimum of ~ 0.32 mag which 
corresponds to the foreground extinction. 

The reddening for the individual star having a photo- 
metric error in the V band oy < 0.1 mag and of a spectral 
type earlier than AO has also been estimated by using 
the reddening free index Q (Johnson & Morgan 1953). 
Assuming the normal reddening slope we can construct 
the reddening-free index Q = (U - B) - 0.72 x (B - V). 
For stars earlier than AO, the value of Q will be less than 
0.0. For main-sequence (MS) stars, the intrinsic (B — V)o 
colour and colour-excess can be obtained from the rela- 
tion (B - V) = 0.332 x Q (Johnson 1966; Hillenbrand et 
al. 1993) and E(B — V) = (B — V) — (B — V)o, respec- 
tively. The individual reddening of stars down to the AO 
spectral type is found to vary in the range of E(B — V)~ 
0.32 - 0.52 mag, implying the presence of a small amount 
of differential reddening. 
4-2.2. Reddening law 

The extinction in star clusters arises due to two distinct 
sources; (i) the general interstellar medium (ISM) in the 
foreground of the cluster, and (ii) the localised dust asso- 
ciated with the cluster. While for the former component a 
value of i? = 3.1 is well accepted (Wegner 1993; Lida et al. 
1995; Winkler 1997), for the intra-cluster extinction the R 
value varies from 2.42 (Tapia et al. 1991) to 4.9 (Pandey 
et al. 2000 and references therein) or even higher depend- 
ing upon the conditions occurring in the region. 
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To study the nature of the extinction law in the IC 1590 
region, we used TCDs as described in Pandey et al. (2000, 
2003). The TCDs of the form of (V - A) vs. (B — V), 
where A is one of the colour bands R, I, J, H, K and L, 
provide an effective method for separating the influence 
of the possible abnormal extinction arising within intra- 
cluster regions having a peculiar distribution of dust sizes 
from that of the normal extinction produced by the diffuse 
interstellar medium (cf. Chini & Wargau 1990, Pandey 
et al. 2000). The TCDs for the nearby reference region 
well away from IC 1590 (see §4.1.2) yield slopes of the 
distributions for (V—I),(V—J),(V—H),(V—K) vs. (B — 
V) as 1.06 ±0.03, 1.97 ±0.05, 2.50 ±0.06, and 2.68 ±0.06 
respectively, manifesting a normal reddening law for the 
foreground interstellar matter (cf. Pandey et al. 2000). 

The (V - K) vs. (B - V) TCDs for the reference and 
cluster regions are shown in Fig. 5. The contamination 
due to field stars is apparent in the cluster region. We 
selected probable field stars having (B — V) > 0.7 visu- 
ally, assuming that stars following the slope of the dis- 
tribution of the reference region are contaminating fore- 
ground stars in the cluster region, and they are shown by 
filled circles. The slopes of the distributions for the prob- 
able cluster members (open circles) , m c i uster are found to 
be 1.24 ±0.04, 2.16 ±0.06, 2.76 ±0.08, 2.90 ±0.09 for the 
(V - I),{V - J), (V - H),(V - A') vs. (B - V) TCDs re- 
spectively. The ratios ^gZy) an< ^ * nc ra ti° of the total- 
to-selective extinction in the cluster region, R c iusten is 
then derived using the procedure given by Pandey et al. 
(2003). Rduster has turned out to be 3.5 ± 0.3. From 
the photometry of bright cluster members in the I and K 
bands and a variable-extinction analysis of ZAMS mem- 
bers, Guetter & Turner (1997) also found a value of R as 
3.44 ±0.07, which is comparable to ours. Several studies 
have already pointed out anomalous reddening laws with 
high R values in the vicinity of star forming regions (see 
e.g. Pandey et al. 2003 and references therein). The 
higher than normal values of R have been attributed to 
the presence of larger dust grains. There is evidence that 
within dark clouds accretion of ice mantles on grains and 
their coagulation due to collision changes the size distri- 
bution towards larger dusts. 

4- 3. Identification of YSOs 

4.3.1. On the basis of (J - H) /(H - A) TCD 

NIR imaging surveys are a powerful tool to detect 
YSOs in star forming regions. The locations of YSOs on 
(J — H)/[H — K) two-colour diagrams (NIR TCDs) arc 
determined to a large extent by their evolutionary state. 
Protostcllar-like objects, CTTSs, weak-line T Tauri stars 
(WTTSs), Herbig Ae/Be stars, and classical Be stars tend 
to occupy different regions on NIR TCDs. 

The NIR TCD using the 2MASS data for all the sources 
lying in the NGC 281 region and having photometric er- 
rors less than 0.1 magnitude is shown in the left panel 
of Fig. 6. All the 2MASS magnitudes and colours have 
been converted into the California Institute of Technology 
(CIT) system. The solid and thick dashed curves rep- 



resent the unreddened MS and giant branch (Bessell & 
Brett 1988) respectively. The dotted line indicates the lo- 
cus of unreddened CTTSs (Meyer et al. 1997). All the 
curves and lines are also in the CIT system. The parallel 
dashed lines are the reddening vectors drawn from the tip 
(spectral type M4) of the giant branch ( "upper reddening 
line"), from the base (spectral type AO) of the MS branch 
( "middle reddening line" ) and from the tip of the intrinsic 
CTTS line ("lower reddening line"). The extinction ratios 
Aj/A v = 0.265, A H /A V = 0.155 and A K /A V = 0.090 have 
been taken from Cohen ct al. (1981). We classified sources 
according to three regions in this diagram (cf. Ojha et 
al. 2004a). 'F' sources are located between the upper and 
middle reddening lines and are considered to be either field 
stars (MS stars, giants) or Class III and Class II sources 
with small NIR-excesses. 'T' sources are located between 
the middle and lower reddening lines. These sources are 
considered to be mostly CTTSs (or Class II objects) with 
large NIR-excesses. There may be an overlap of Herbig 
Ae/Be stars in the 'T' region (Hillenbrand et al. 1992). 
'P' sources are those located in the region redward of the 
lower reddening line and are most likely Class I objects 
(protostar-like objects; Ojha et al. 2004a). It is worth- 
while to mention also that Robitaille et al. (2006) have 
shown that there is a significant overlap between proto- 
stars and CTTSs. The NIR TCD of the NGC 281 region 
(left panel of Fig. 6) indicates that a significant number 
of sources show (H — K) excess and these are shown by 
open triangles. The sources having X-ray emission and 
Ha emission are shown by circles (open and filled) and 
star symbols respectively. A comparison of the TCD of 
the NGC 281 region with the NIR TCD of nearby ref- 
erence region (right panel of Fig. 6) indicates that the 
sources in the NGC 281 region having X-ray emission and 
lying in the 'F' region above the extension of the intrin- 
sic CTTS locus as well as sources having (J — H) > 0.6 
mag and lying to the left of the first (left-most) reddening 
vector (shown by filled circles) could be WTTSs/Class III 
sources. Here it is worthwhile to mention that some of 
the X-ray sources classified as WTTSs/ Class III sources, 
lying near the middle reddening vector could be CTTSs/ 
Class II sources. The CTTSs and WTTSs identified in 
this section are listed in Table 7. 
4.3.2. On the basis of MIR data 

The NGC 281 region also has MIR observations through 
the Spitzer Space Telescope towards the south-west direc- 
tion of the cluster. Since young stars inside cloud clumps 
are often deeply embedded, these MIR observations can 
provide a deeper insight into the embedded YSOs. YSOs 
occupy distinct regions in the IRAC colour plane accord- 
ing to their nature; this makes MIR TCDs a very useful 
tool for the classification of YSOs. Whitney et al. (2003) 
and Allen et al. (2004) presented independent model pre- 
dictions for IRAC colours of various classes of YSOs. Fig. 
7 (left) presents a [5.8]- [8.0] versus [3.6]- [4.5] TCD for the 
observed sources. The sources within the box represent 
the location of Class II objects (Allen et al. 2004; Megeath 
et al. 2004). The sources located around [5.8]-[8.0] = 
and [3.6]-[4.5] = are foreground/background stars as 
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well as Class III objects. Sources with [3.6]- [4.5] > 0.8 
and/or [5.8]-[8.0] > 1.1 have colours similar to those de- 
rived from models of protostcllar objects with in-falling 
dusty envelopes (Allen ct al. 2004). These are Class 0/1 
sources. Encircled objects represent sources with X-ray 
emission. A majority of the Class II objects have X-ray 
emission, whereas none of the Class 0/1 sources show x- 
ray emission. It is found that four of the probable Class 
0/1 sources identified on the basis of MIR data lie in the 
unexposed area of the detector (ACIS-I) of the Chandra 
telescope, whereas one lies near the edge of detector ACIS- 

I of the Chandra telescope. 

The detection of Class 0/1 and Class II sources in all 
the four IRAC bands is limited mainly by the lower sensi- 
tivity of the 5.0 and 8.0 fxm channels. Fig. 7 (right panel) 
shows the IRAC colour-magnitude diagram (CMD) for 
stars detected in the 3.6 and 4.5 fim bands. Encircled 
objects represent sources with X-ray emission. Stars hav- 
ing 0.35 < [3.6] - [4.5] < 0.80 mag could be probable Class 

II stars, whereas stars having [3.6]- [4.5] > 0.80 could be 
Class 0/1 sources. As can be seen, a majority of the Class 
II sources are X-ray emitting stars, whereas X-ray emis- 
sion is mostly absent in probable Class 0/1 sources. One of 
the critical astrophysical questions is whether X-ray emis- 
sion is present in Class /I sources at the very onset of star 
formation when collimatcd outflows begin (Gctman ct. al 
2007). A few studies report detection of X-ray from Class 
0/1 protostars, whereas some studies reported that many 
bona-fide protostars arc not detected in X-ray images (cf. 
Tsuboi et al. 2001, Hamaguchi et al. 2005, Getman et al. 
2007). The non-detection of X-ray emission in Class 0/1 
sources is usually attributed to heavy obscuration instead 
of the intrinsic absence of X-ray emission in protostars. 
The sources having colours 0.35 < [3.6] — [4.5] < 0.80 mag 
and X-ray emission are considered as Class II sources and 
these are also listed in Table 7. 

To further elucidate the nature of the PMS sources, we 
derived SEDs for 35 sources (cf. Table 8) using the optical, 
NIR and MIR photometry. In Fig. 8, we show a sample of 
three SEDs for three different classes. To classify the evo- 
lutionary stage of YSOs using the SEDs, we adopted the 
classification scheme of Lada ct al. (2006), which defines 
the spectral class index a = (flog A(7\)/<flog(A). We com- 
puted the spectral class index olk— 8/mu which is the slope 
of the linear fit to the fluxes at the K s and IRAC 8 fj,m 
bands. Objects with axs^m > +0.3, +0.3 > aK-8fim > 
-0.3, -0.3 > OLK-i^m > —1.8 and -1.8 > a K ^ m are 
considered as Class I, Flat, Class II and Class III sources 
respectively. The a indices obtained from the SEDs, in 
general, confirm the classification obtained from the MIR 
TCD (cf. Fig. 7 (left)). However, the ax-s^m indices 
of two sources (Ih and Ij; Class I on the basis of MIR 
TCD) reveal that these must be Class II sources, whereas 
two seemingly Class II sources on the basis of the MIR 
TCD (Ha and lie) appear to be Class I sources on the ba- 
sis of the SEDs. The MIR TCD classifies the two sources, 
namely Illb and Illh as Class III sources, whereas oiK-s^m 
indices classify them as Class II objects. 

Table 7 provides a complete list of YSOs identified in 



the present study on the basis of Ha emission, X-ray emis- 
sion, NIR and MIR observations. The table contains 12 
Ha emission, 134 X-ray emission, 87 NIR excess (CTTSs) 
source and 118 WTTSs. The MIR data yield 25 and 61 
Class I and Class II sources, respectively. The J/ (J — 77) 
CMD (Fig. 9) reveals that the identified YSOs are prob- 
ably PMS stars of age < 1 Myr. Majority of these stars 
have masses between 0.5-3.5 M . 

4-4- Optical colour-magnitude diagram 

The V I (V — I) CMD for stars lying within the cluster 
region is shown in the left-hand panel of Fig. 10. A more 
or less well defined broad MS presumably due to the vari- 
able reddening in the cluster region can be noticed down 
to ~15 mag. The distribution of stars fainter than V ~ 
15-16 mag deviates towards the red side of the MS in- 
dicating the presence of PMS stars in the cluster region. 
Contamination due to a field star population is also evi- 
dent in the CMD. To study the LF/MF of the cluster, it 
is necessary to remove field star contamination from the 
sample of stars in the cluster region because PMS mem- 
ber stars and dwarf foreground stars both occupy similar 
positions above the ZAMS in the CMD. In the absence of 
proper motion data, we used a statistical method to esti- 
mate the number of probable member stars in the cluster 
region. We again utilize the reference region towards the 
north-west (cf. §4.1.2); it has the same area as that of 
the cluster region. The middle panel of Fig. 10 shows its 
V/(V-I) CMD. 

To remove contamination due to field stars, we statisti- 
cally subtracted their contribution from the CMD of the 
cluster region using the following procedure. The CMDs 
of the cluster as well as of the reference region were divided 
into grids of AV = 1 mag by A(V — I) = 0.4 mag. The 
number of stars in each grid of the CMDs were counted. 
After applying the completeness corrections using the CF 
(cf. Table 2) to both the data samples, the probable 
number of cluster members in each grid were estimated 
by subtracting the corrected reference star counts from 
the corrected counts in the cluster region. The estimated 
numbers of contaminating field stars were removed from 
the cluster CMD in the following manner. For a ran- 
domly selected star in the CMD of the reference region, 
the nearest star in the cluster CMD within V ± 0.25 and 
(V — 7) ± 0.125 of the field was removed. Although the 
statistically cleaned V/ (V — I) CMD of the cluster region 
shown in Fig. 10c clearly shows the presence of PMS stars 
in the cluster, however the contamination due to field stars 
at V > 17 mag and (V — 7) ~ 1.2 mag can still be seen. 
This field population could be due to the background pop- 
ulation as discussed by Pandey et al. (2006). 

The available distance estimates of NGC 281 in the lit- 
erature varies from 2.0 kpc to 3.7 kpc (cf. Sato et al. 
2008). For further analysis we adopt the VLBI trigono- 
metric distance of the maser source of 2.81 ± 0.24 kpc 
(Sato et al. 2008). Fig. 11 (left panel) shows statistically 
cleaned dereddened Vo/(V — I)q CMD where stars having 
spectral type earlier than AO were individually dereddened 
(cf. §4.2.1), whereas the mean reddening of the nearby 
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region, estimated from the available individual reddening 
values in that region, was used for other stars. We have 
also plotted the ZAMS by Marigo ct al. (2008) and the 
PMS isochrones by Siess et al. (2000) using the distance 
of 2.81 ±0.24 kpc. The evolutionary tracks by Siess et al. 
(2000) for various masses have also been plotted which 
reveal that majority of the YSOs have masses between 
0.5-3.5 M©. Fig. 11 (left panel) indicates an age spread 
for the PMS population. To check its reality, we plot- 
ted V /(V - J) CMD (assuming a mean E(B - V)=0A 
mag) for the Ha emission stars, NIR-excess stars (proba- 
ble CTTSs) and X-ray stars (probable WTTSs) (cf. §4.3) 
in Fig. 11 (right panel). This also indicates an age spread 
of about 1-5 Myr for these probable PMS stars, support- 
ing the reality of that in Fig. 11 (left panel). 

The age and mass of each YSO were estimated by com- 
paring its location with the isochrones. Here we would like 
to point out that the estimation of the ages and masses 
of the PMS stars by comparing their positions in the 
CMDs with the theoretical isochrones is prone to random 
as well as systematic errors (see e.g. Hillenbrand 2005, 
Hillenbrand 2008, Chauhan et al. 2009, 2011). Chauhan 
et al. (2009) and Barentsen at al. (2011) have studied 
the effect of random errors in the age estimation of PMS 
stars. Barentsen at al. (2011) found that uncertainty 
in the extinction estimation could play a significant role. 
In the case of NGC 281 region the variable extinction is 
small (~ 0.2 mag), hence it should not contribute signif- 
icantly in the errors. The effect of random error due to 
photometric error and reddening estimation in determina- 
tion of ages and masses was estimated by propagating the 
random errors to their observed estimation by assuming 
normal error distribution and using the Monte-Carlo sim- 
ulations. The estimated ages and their error are given in 
Table 9. The systematic errors could be due to the use of 
different PMS evolutionary models and the error in dis- 
tance estimation etc. Barcntson et al. (2011) mentioned 
that the ages may be wrong by a factor of two due to the 
systematic errors in the model. The presence of binaries 
may be the another source of error in the age determi- 
nation. Binarity will brighten the star, consequently the 
CMD will yield a lower age estimate. In the case of an 
equal mass binary we expect an error of ~ 50 - 60% in the 
PMS age estimation. However, it is difficult to estimate 
the influence of binaries/ variables on mean age estimation 
as the fraction of binaries/variables is not known. In the 
study of TTSs in the HII region IC 1396, Barentson et al. 
(2011) presumed that the number of binaries in their sam- 
ple of TTSs could be very low as close binary lose their 
disc significantly faster than single stars (cf. Bouwman ct 
al. 2006). 

The age distribution of YSOs shown in the Fig. 12 indi- 
cates a significant scatter. The inset of Fig. 12 shows the 
distribution of random errors. A comparison manifests 
that the age distribution of YSOs shows a significantly 
larger scatter than that could occur due to random errors. 
Burningham et al. (2005) have investigated the effect of 
photometric variability in the apparent age spreads ob- 
served in the CMDs of OB associations. They found that 



the combination of binarity, photometric uncertainty and 
variability could not explain the observed age spread in 
the CMDs of OB associations. If the effect of unresolved 
binaries is not significant, we presume that the main rea- 
son of spread in the distribution of YSOs in Fig. 11 could 
be due to the spread in the ages of YSOs. 

The ages of young clusters are usually derived on the 
dereddened CMDs by comparing the earliest members 
to post-main-sequence evolutionary tracks if significant 
evolution has occurred and/or the low mass contracting 
population to the PMS isochrones. Since the most mas- 
sive member of the cluster IC 1590 is an 06.5 MS star 
(Walborn 1973; Abt 1986; Guetter & Turner 1997), its 
maximum age should be of the order of the MS life time 
of this star, i.e., ~4.4 Myr (Meynet et al. 1994). Based 
also on the PMS stars, Guetter & Turner (1997) derived 
a maximum age of ~ 3.5 Myr for IC 1590. 

We consider the points lying above the 5 Myr isochronc 
in Fig. 11 (left panel) as representing the statistics of the 
PMS stars in the cluster region. Here, we would like to 
remind the readers that the filled circles in Fig. 11 may 
not represent the actual members of the cluster. However, 
they should represent the statistics of PMS stars in the 
cluster region and this statistics is used to study its MF 
only. 

4-. 5. Initial Mass Function and K-band luminosity func- 
tion 

The distribution of stellar masses that form in a star- 
formation event in a given volume of space is called initial 
mass function (IMF). Young clusters are important tools 
to study the IMF since they are too young to lose a sig- 
nificant number of members cither by dynamical or by 
stellar evolution. The MF is often expressed by a power 
law, Y(logm) cx m T and the slope of the MF is given as: 

r = dlogiV(logm)/dlogm 

where Y(logm) is the number of stars per unit logarith- 
mic mass interval. The classical value derived by Salpeter 
(1955) is T = -1.35. 

With the help of the statistically cleaned CMD, shown 
in Fig. 11 (left panel), we can derive the MF using the 
theoretical evolutionary models. Since the age of the mas- 
sive cluster members is thought to be ~ 2 - 4 Myr, the 
stars having V < 13 mag (Vq < 12 mag; M > 4M Q ) are 
considered to be still on the MS. For these stars, the LF 
was converted to a MF using the theoretical models by 
Marigo al. (2008) (cf. Pandey et al. 2001, 2005). The 
data for the three brightest stars, which were saturated 
in the present photometry, have been taken from Guetter 
& Turner (1997). The MF for the PMS stars have been 
obtained by counting the number of stars in various mass 
bins (shown as evolutionary tracks) having age < 5 Myr 
in Fig. 11 (left panel). The resulting MF of the cluster 
is plotted in Fig. 13. Since data incompleteness plays an 
important role in the estimation of the IMF, we restrict 
our analysis only to sources having V < 16.5 mag. The 
present data have a completeness of ^95% at V =16.5 
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mag (cf. Table. 2). The slope (T) of the MF in the 
mass range 2 < M/Mq < 54 comes out to be —1.11 ±0.15, 
which seems to be slightly shallower than the Salpetcr 
(1955) value (-1.35). Using various combinations of the 
maximum expected errors in E(B — V) and the distance, 
we found that the slope of the MF can vary in the range 
of -1.08 ±0.15 to -1.16 ±0.15. Keeping the errors in the 
estimation of T in mind, it is difficult to decide whether 
the present slope T' is different from that of the Salpcter 
value. Guetter & Turner (1997) have reported the slope 
of the MF (r = -1.00 ±0.21) for IC 1590 which, within 
error, is comparable to the value obtained in the present 
work. 

The K-band luminosity function (KLF) is a powerful 
tool to investigate the IMF of young embedded clusters; 
therefore during the last decade several studies focused 
on the determination of the KLFs of young open clusters 
(e.g. Lada & Lada 2003, Ojha et al. 2004b, Sanchawala 
ct al. 2007). In order to obtain the KLF of IC 1590, we 
again have to examine the effects of incompleteness and 
field star contamination in our data. The completeness 
of the data is estimated using the ADDSTAR routine of 
DAOPHOT as described in §2.1.1. To take into account 
the foreground/background field star contamination, we 
used the Besangon Galactic model of stellar population 
synthesis (Robin et al. 2003) and predicted the star counts 
in both the cluster region and in the direction of the refer- 
ence field. We checked the validity of the simulated model 
by comparing the model KLF with that of the reference 
field and found that the two KLFs match rather well (Fig. 
14a) . An advantage of using the model is that we can sep- 
arate the foregrounds (d < 2.8 kpc) and the background 
(d > 2.8 kpc) field stars. The foreground extinction to- 
wards the cluster region is found to be Ay ~ 1.0 mag. The 
model simulations with d < 2.8 kpc and Ay — 1.0 give 
the foreground contamination, and that with d > 2.8 kpc 
and Ay = 1.7 mag the background population. We thus 
determined the fraction of the contaminating stars (fore- 
ground+background) over the total model counts. This 
fraction was used to scale the nearby reference region and 
subsequently the modified star counts of the reference re- 
gion were subtracted from the KLF of the cluster to ob- 
tain the final corrected KLF. This KLF is expressed by 
the following power-law: 

dN(K) i naif 
dK 01 1U 

where d ^|^ is the number of stars per 0.5 magnitude bin 
and a is the slope of the power law. Fig. 14b shows 
the KLF for the cluster region. This indicates a slope 
of a = 0.37 ±0.07 which is similar to the average slopes 
(a ~ 0.4) for young clusters (Lada et al. 1991; Lada & 
Lada 1995; Lada & Lada 2003) but higher than the values 
(0.27 - 0.31) obtained for Be 59 (Pandey et al. 2008) and 
Stock 8 (Jose et al. 2008). 

5. Star formation scenario 

The star forming region NGC 281 has been attract- 
ing attention of the star forming community. It contains 



a cluster IC 1590 with a Trapezium like system of O- 
type stars at the center. The ionized hydrogen seems to 
be associated with two CO molecular clumps (east and 
west) which have been mapped in 12 CO and 13 CO by 
Elmegreen & Lada (1978), Leisawitz, Bash & Thaddeus 
(1989), Henning et al. (1994), Megeath & Wilson (1997) 
and Lee & Jung (2003). The western CO clump called 
NGC 281 West is somewhat more massive and compact 
than the elongated eastern clump (NGC 281 East) (Lee 
& Jung 2003). 

The ionizing source HD 5005 lies to the north- 
east/northwest of NGC 281 West/NGC 281 East. The 
differential extinction towards the central cluster is ~0.2 
mag (cf. §4.2.1) indicating that the central cluster con- 
tains only gas and dust of low-density. A similar trend has 
been noticed in many clusters associated with HII regions 
(e.g. 30 Dor, Brandl et al. 1996; NGC 3603, Pandey et 
al. 2000 and NGC 1893, Sharma et al. 2007). A reason- 
able explanation for this lack of a dense medium in the 
central region may be the effects of the massive star(s) at 
the cluster center. 

As indicated in previous studies (cf. Lee & Jung 2003 
and references therein), the western molecular clump is 
interacting with the ionized gas. On the basis of kinetic 
evidence, Elmegreen & Moran (1979) suggested the pas- 
sage of a shock through NGC 281 West. An H2O mascr 
was found to be coincident with the peak of the cloud in- 
dicating the ongoing star formation (Elmegreen & Lada 
1978). Carpenter et al. (1993) and Megeath (1994) de- 
tected a cluster of low mass stars associated with NGC 
281 West. Elmegreen & Lada (1978) suggested that this 
region is a site of triggered star formation through the 
"collect and collapse" mode. However, Megeath & Wilson 
(1997) claimed that numerical models of imploding spher- 
ical clumps can approximately reproduce the kinematic 
features observed in NGC 281 West, and suggested that 
"radiation-driven implosion" (RDI) is a more plausible 
and attractive model. Between NGC 281 East and the 
central cluster a few bright-rimmed clouds (BRCs) or 
cometary globules and an IRAS source are located. Some 
of the Ha stars can be seen around the tip of the BRCs 
just as observed in several well-known BRCs (see Ogura 
et al. 2002). Both the CO clumps show the presence of 
YSOs (IR-excess and Ha stars) around their respective 
centers. 

Whatever the star formation scenario, the ioniza- 
tion/shock fronts caused by high- mass stars of the first 
generation appears to have initiated the formation of a 
new generation of stars at the edge of the molecular 
clumps. The distribution of YSOs and morphological de- 
tails of the environment around the cluster can be used to 
infer the star formation history of the NGC 281 region in 
detail. To know the distribution of YSOs in the region, we 
have included all the detected YSOs without considering 
their photometric errors to improve the sample. 

IRAS maps can be used to study the distribution of 
dust and unidentified infrared band (UIB) carriers. Fig. 
15 shows the IRAS intensity maps for the NGC 281 re- 
gion at 12 fim (top left), 25 fim (top right), 60 /^m (bottom 
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left) and 100 fim (bottom right). The global features of 
these maps are quite similar to that of the molecular gas. 
All the contours show peaks around both the eastern and 
western clumps indicating the coexistence of warm dust 
(IRAS 12 and 25 /mi), cold dust (IRAS 60 and 100 /xm) 
and molecular gas. NGC 281 East lacks radio continuum 
emission and its extended distribution of the 12 /mi emis- 
sion towards the south resembles that of the molecular 
material. As pointed out by Leisawitz, Bash & Thaddeus 
(1989), in the case of eastern clump, the peak of ionized 
gas (as seen from the ionized source) is followed by the 
peaks of IR emission and CO emission from the molecular 
cloud, respectively. The coincidence of the peak of the 
IRAS and CO emission, the location of the IRAS point 
source and the distribution of YSOs indicate ongoing star 
formation activity in these clumps. 

Fig. 16 shows the map of the 12 CO emission taken 
from Henning et al. (1994) and the 1.4 GHz radio emis- 
sion from NVSS along with the spatial distribution of all 
the detected YSOs overlaid on the DSS-2 R band image. 
The center of the cluster IC 1590 is marked by 'C and 
the location of the ionization source is shown by a square. 
A well aligned distribution of the detected YSOs from the 
vicinity of the ionization source to the direction of NGC 
281 West can be noticed. This spatial distribution of the 
YSOs resembles with that in the case of NGC 1893 and 
BRC 14, where a similar distribution of NIR-excess stars 
can be noticed from the ionization source to the direction 
of the cometary globules/BRC (see figure 22 of Sharma 
et al. 2007 and figure A3 of Chauhan et al. 2009). These 
alignments in NGC 1893 and BRC 14 were attributed to 
triggered star formation due to a series of RDI process. 
In both the cases YSOs located away from the ionization 
sources are found to be younger. Sicilia-Aguilar ct al. 
(2004) have also shown that in the case of the Tr 37/IC 
1396 globule region, CTTSs are aligned from the iosizing 
source towards the direction of the globule, and that most 
of the younger (^1 Myr) members appear to lie near or 
within the globule. They mentioned that it can be in- 
dicative of triggered star formation. Fig. 17 (left panel) 
shows the age distribution of the YSOs as a function of 
the radial distance from HD 5005, the ionization source 
of NGC 281, which shows that the YSOs in the cluster 
region (r < 8') have ages < 5 Myr whereas those lying 
outside the cluster region are relatively younger and have 
ages < 2 Myr. Fig. 17 (right panel) shows the radial 
variation of the NIR-excess A(H — K), which is defined 
as the horizontal displacement from the middle redden- 
ing vector at the boundary of the 'F' and 'T' regions (see 
Fig. 6). To quantify the radial variation of the age and 
NIR-excess, we used only those stars which have the error 
< 0.1 mag. The distribution of A(H — K) also suggests 
that the sources lying outside the boundary of the cluster 
(r ~ 8') have relatively larger NIR excess in comparison to 
those located within the cluster region. However we ad- 
mit that the above statements are not conclusive in view 
of the fact that the scatters are large and the differences 
are subtle. A deeper optical, NIR and MIR observations 
are needed to have conclusive star formation scenario in 



the region. A similar trend has been reported in the case 
of a few BRCs by Chauhan et al. (2009) as well as in a 
recent study on the IC 1396 region by Barcntscn et al. 
(2011). 

The near-IR excess in the case of CTTSs suggest the 
presence of dusty optically thick discs (Haisch et al. 2001a, 
Sicilia-Aguilar et al. 2006). Sicilia-Aguilar et al. (2006) 
found similarities between the decrease in IR excesses and 
the decrease in accretion rates, and concluded that gas 
evolution seems to occur somehow parallel to the evolu- 
tion of the dust grains and the structure of the disc. In 
a recent study, Sicilia-Aguilar ct al. (2009) further con- 
firmed that IR-excess decreases with age. Hence, rela- 
tively large NIR excess outside the cluster region (cf. Fig. 
17, right panel) could be indicative of relatively younger 
population outside the cluster region. 

To study the evolutionary stages of the NGC 281 re- 
gion, we have divided it into three sub-regions as shown 
in Fig. 16, namely Region 1 (IC 1590), Region 2 (NGC 
281 West) and Region 3 (NGC 281 East). Figs 18 and 19 
show the V/(V-I) CMDs and NIR TCD, respectively, for 
the three regions. Here also we have used only sources of 
good quality, i.e., sources having error less than 0.1 mag. 
X-ray data are not available for Region 3. The CMDs indi- 
cate that the age of the YSOs in Region 1 ranges between 
< 1 - 5 Myr, whereas a majority of the YSOs in Region 

2 have ages < 1 Myr. The YSOs associated with Region 

3 indicate an age of ~ 1-2 Myr for the region. The NIR 
TCDs (Fig. 19) indicate that the extinction in Regions 2 
and 3 is relatively higher in comparison to Region 1. 

A comparison of the statistics of the YSOs in these three 
regions can also give a clue about the evolutionary stages 
of the regions. Table 10 shows the statistics of the prob- 
able CTTSs associated with the three sub-regions. We 
have given statistics for the sources having error < 0.1 
mag as well as for all the detected YSOs without consid- 
ering their errors. The total number of stars detected in 
the 2MASS Catalogue and the number of field stars ex- 
pected in each region, estimated from the nearby reference 
region (cf. §4.1.2), are also given in Table 10. The per- 
centage of Ha stars and NIR excess stars in each region is 
estimated after substracting the contribution of expected 
field stars (89 stars, column 2 of Table 10) from the total 
stars (also given in column 2 of Table 10) in the subre- 
gions. The fraction of detected CTTSs (i.e. NIR-excess 
stars and Ha stars), for both the samples, is significantly 
higher in Regions 2 and 3 as compared to the cluster re- 
gion (Region 1). This is further supported by a compari- 
son of the CTTSs fractions, }ctts = Nctts/ (Notts + 
Nwtts ) in Region 1 (5-12 %) and Region 2 (25-45 %) 
(cf. Table 11). The sources flagged as 1 and 2 in column 
14 of Tabic 7 arc classified as CTTSs and WTTSs respec- 
tively. The WTTSs are identified using the NIR TCD 
and X-ray data (cf. §4.3.1). The X-ray data is not avail- 
able for the Region 3. Here it is important to mention 
that as mentioned in §4.3.1 some of the CTTSs having 
less or negligible NIR excess might have been classified as 
WTTSs. If this is true, the fcTTS will have further higher 
value. Haisch et al. (2001b) have found the disc evolu- 
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tion fraction in the sense that the fraction of stars having 
a disc decreases with age. Armitage et al. (2003) have 
also found that, in Taurus-Auriga T-association, /ctts, 
decreases with stellar age. A comparison of the disc frac- 
tions of Regions 1 and 2 with those given by Haisch et 
al. (2001b) and Armitage et al. (2003) suggests ages of ~ 
4-5 Myr and ~ 1-3 Myr for Regions 1 and 2 respectively, 
which is in fair agreement with the results obtained from 
the CMDs. Thus the above discussions suggest the prop- 
agation of (triggered) star formation in Regions 2 and 3. 

Megeath & Wilson (1997) pointed out the presence of 
two sub-clusters (northern and southern) in NGC 281 
West. They concluded that the northern and southern 
sub-clusters are indeed physically separate sub-clusters re- 
sulting from distinct star formation events. Based on the 
location of the sub-clusters they concluded that the north- 
ern sub-cluster appears to be associated with the NW and 
NE clumps, whereas the southern sub-cluster seems to be 
associated with S clump. The spatial distribution of the 
Class 0/1 and Class II sources detected by Spitzer obser- 
vations overlaid on the 2MASS image is shown in Fig. 20. 
Certainly the YSOs are found to make two sub-clusters. 
However, a careful look of this figure manifests a third sub- 
clustering also. The isodensity contours shown in Fig. 2 
also suggest the presence of three sub-clusters. We have 
visually marked the boundaries of these sub-clusters, des- 
ignated as 'a', 'b' and 'c'. Some of the detected YSOs 
are located outside the boundaries of these sub-clusters. 
Megeath & Wilson (1997) also found that northern sub- 
cluster has higher extinction as compared to the south- 
ern sub-cluster. A comparison of the morphology of the 
molecular cloud as observed in C 18 (2-1) (Megeath & 
Wilson 1997) and the distribution of YSOs in the 'a' sub- 
cluster region indicates that 'a' sub-cluster seems to be 
associated with the NW and NE clumps as referred by 
Megeath & Wilson (1997), whereas 'c' sub-cluster seems 
to be associated with the clump 'S'. Megeath & Wilson 
(1997) found that the clump 'S' (-34 to -32 kms" 1 ) is 
kincmatically and spatially distinct from the NW and NE 
clumps (-32 to -29 kms -1 ). The velocity distribution sug- 
gests that clump 'S' (i.e. sub-cluster 'c') is relatively near 
to the observer in comparison to the NW clump 'a'. If 
we assume that the sub-clusters are not associated with 
each other, the spatial distribution of the YSOs in the 
sub-cluster 'a' reveals that the Class II sources are rel- 
ativley near to the ionizing source in comparison to the 
Class 0/1 sources. We do not find this trend in the sub- 
cluster 'b'. The distribution of YSOs in sub-cluster 'c' 
also show a similar trend however the statistics is poor. 
The distribution of the YSOs detected by using the IRAC 
data in the case of a few BRCs also shows that the Class 
0/1 sources are found to be located more away from the 
ionizing sources as compared to Class II sources. This 
seems to indicate the propagation of star formation in the 
NGC 281 West region. 



6. PMS sources: MIR and X-ray observations 

As discussed in §4.3.2 that none of the Class 0/1 sources 
show X-ray emission. The high value of Ay in the region 
could be a possible reason for non-detection of X-ray in 
the Class 0/1 sources. X-ray median energy (MedE) is a 
reliable indicator of absorbing column density through the 
empirical relation, log A H = 21.22 + 0.44 (MedE) cm~ 2 
(Fcigelson et al. 2005). Getman et al. (2007) have demon- 
strated a relationship between MedE and the MIR colour 
[3.6]-[4.5] for obscured PMS stars. Fig. 21 shows a MedE 
vs. [3.6]-[4.5] diagram. The dashed lines represent the 
boundaries of the distribution obtained in IC 1396 N by 
Getman et al. (2007, see their figure 8). Using the Spitzer 
photometry and assuming that the relationship between 
MedE and MIR colours given by Getman et al. (2007) is 
valid for the NGC 281 region also, we estimated the MedE 
for the Class II /III sources detected in the present study 
as < 4 keV. Using the above relation, the absorbing col- 
umn densities towards these sources have been estimated 
as Njj > 10 23 cm~ 2 , which yields Ay > 50 mag according 
to the standard gas-to-dust ratio by Ryter (1996). Here it 
is interesting to mention that all the six Class /I sources 
detected in the IC 1396 N region (distance ^0.75 kpc) by 
Getman et al. (2007) having Ay > 50 mag show X-ray ac- 
tivity. Non-detection of X-ray emission in probable Class 
0/1 sources in the NGC 281 region may be due to the de- 
tection limit on account of the larger distance. Deeper 
exposures are required to reach the conclusions. 

7. Summary and conclusions 

In this paper, we present a multiwavelength study of 
the NGC 281 region using deep wide-field optical U BVI C 
data, slitless spectroscopy along with the archival data 
from the surveys such as Chandra, Spitzer, 2MASS, IRAS 
and NVSS. We made an attempt to construct a scenario 
of the global star formation in the NGC 281 complex by 
taking into account the effects of massive stars on low- 
mass star formation. The main results from this study 
are as follow: 

• The morphology of the central cluster IC 1590 is 
found to be elongated. The extent of the cluster is 
6.5 pc and the core radius 1.6 pc. The maximum 
age of the ionizing source of the region could be ~ 
4 Myr. The minimum reddening E(B — V) m i„ to- 
wards the cluster is estimated to be 0.32 mag and the 
cluster shows a small amount of differential redden- 
ing (E(B — V) ~ 0.2 mag). The distribution of the 
YSOs selected on the basis of NIR-excess, Ha emis- 
sion and X-ray emission show a spread in the CMD. 
The age distribution of YSOs indicates a non-coeval 
star formation in and around the cluster. 

• A slitless spectroscopic survey of the NGC 281 re- 
gion identifies 12 Ha emission stars. Some of these 
stars are located near the globulcs/BRCs and show 
the properties of intermediate-low mass PMS stars. 
On the basis of NIR excess we identified 87 CTTSs. 
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Using a Chandra archival dataset and NIR colour- 
colour diagram we also identified 118 WTTSs. A 
majority of the identified YSOs (IR-excess, X-ray 
and Ha stars) are low mass PMS stars having ages 
< 1- - 2 Myr and masses 0.5-3.5 M Q . 

• The slope (T) of the MF for the central cluster IC 
1590 in the mass range 2 < M/M Q < 54 is found to 
be — 1.11 ±0.15. The slope of the K-band luminosity 
function (0.37 ±0.07) is similar to the average value 
(~0.4) reported for young clusters in the literature 
(i.e. Lada ct al. 1991; Lada & Lada 1995; Lada & 
Lada 2003). 

• The distribution of gas and dust obtained from 
the IRAS, CO and radio continuum maps indicates 
clumpy structures around the central cluster. The 
radial distribution, ages and NIR-excesses A(H — K) 
of the YSOs as well as the fraction of CTTSs, sug- 
gest triggered star formation around the cluster. 
However we would like to caution the readers that 
the above statement is not conclusive in view of 
the scatter in the data. Deeper optical, NIR and 
MIR observations are needed to have a conclusive 
view for the star formation scenario in the region. 
The Class 0/1 and Class II sources detected by us- 
ing the Spitzer MIR observations indicate that a 
majority of the Class II sources are X-ray emitting 
stars, whereas X-ray emission is absent in Class 0/1 
sources. The spatial distribution of Class 0/1 and 
Class II sources reveals the presence of three sub- 
clusters in the NGC 281 West region. The distribu- 
tion of the Class 0/1 and Class II sources in the 'a' 
sub-cluster indicates that the Class II sources tend 
to be located relatively near to the ionizing source. 
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Table 1. Log of observations. 



Date oi observation/ _b liter 


Fxp. (sec) x No. of frames 




Kiso Schmidt telescope, Japan 


21 November 2004 




B 


60 x 6,20 x 6 


V 


60 x 6,10 x 6 


Ic 


60 x 6,10 x 6 


27 November 2005 




U 


180 x 6,60 x 2 


Ic 


10 x 2 




Sampurnanand telescope, ARIES 


07 January 2005 




j j 
u 


oOO X 6, IzU X 1,60 X 1 


D 


120 x 3,30 x 3 


V 


120 x 3,30 x 3 


Ic 


60 x 4,10 x 3 




Himalayan Chandra Telescope, IIA 


10 October 2005 




Slitless spectra 


(420 x 3) x 4 


Direct Frames 


(60 x 3) x 4 


16 August 2006 




Slitless spectra 


(300 x 2) x 4 


Direct Frames 


(60 x 1) x 4 



Table 2. Completeness Factor (CF) of the optical photometric data in the cluster and field regions. 



V range 


IC 1590 




Field region 


(mag) 


r < 2' 2' < r < 5' 


r < 5' 


9.5-10.5 


1.00 


1.00 


1.00 


10.5-11.5 


1.00 


1.00 


1.00 


11.5-12.5 


1.00 


1.00 


1.00 


12.5-13.5 


1.00 


1.00 


1.00 


13.5-14.5 


1.00 


1.00 


1.00 


14.5-15.5 


1.00 


0.96 


0.98 


15.5-16.5 


0.93 


0.97 


0.95 


16.5-17.5 


0.74 


0.94 


0.96 


17.5-18.5 


0.62 


0.90 


0.90 


18.5-19.5 


0.32 


0.40 


0.45 
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Table 3. Comparison of the present photometry with the available photometry in the literature. The differ- 
ence A (literature-present data) is in magnitude. Mean and <r are based on N stars in a V magnitude bin. 



V range 


A(V) 




A(B-V) 




A(U-B) 






(Mean ± a) 


(N) 


(Mean ± a) 


(N) 


( Mean±a) 


(N) 


Gucttcr & Turner (1997, ccd) 














< 12 


-0.004±0.019 


5 


-0.022 ±0.016 


5 


-0.046 ±0.046 


5 


12-13 


0.001 ±0.031 


5 


-0.019 ±0.010 


5 


-0.038 ±0.049 


5 


13-14 


-0.011 ±0.022 


18 


0.002 ±0.018 


18 


-0.040 ±0.075 


16 


14-15 


0.022 ±0.042 


16 


-0.014 ±0.022 


16 


0.010 ±0.086 


14 


15-16 


0.011 ±0.040 


40 


0.011 ±0.035 


40 


0.051 ±0.102 


21 


16-17 


0.021 ±0.053 


52 


-0.026 ±0.042 


50 






17-18 


0.036 ±0.062 


16 


-0.002 ±0.077 


16 






Gucttcr & Turner (1997, pe) 














< 12 


-0.001 ±0.014 


5 


-0.030 ±0.013 


5 


-0.054 ±0.038 


5 


12-13 


-0.018±0.011 


4 


-0.016±0.014 


4 


-0.059 ±0.040 


4 


13-14 


0.015±0.011 


4 


-0.013 ±0.014 


4 


-0.083 ±0.045 


4 



ccd: charged coupled device data 
pe: photo-electric data 



Table 4. The optical, 2MASS and IRAC counterparts of the X-ray sources searched within a match radius of 1 arc- 
sec. The radial distance is from the cluster center. The complete table is available in the electronic form only. 

Radial Distance c*(20oo) ^(2000) V {V ~ -0 J H K s 3.6/im 4.5/im 

(') ( h : m : s ) (°:':"' (mag) (mag) (mag) (mag) (mag) (mag) (mag) 

O07 00:52:39.23 +56:37:49.1 L381 1377TT 131)35 12J43 

0.17 00:52:38.84 +56:37:37.0 14.391 1.029 12.707 12.229 12.115 

0.41 00:52:38.64 +56:38:09.2 19.155 1.837 16.015 15.501 15.047 

0.44 00:52:42.48 +56:37:55.9 - - 16.060 15.143 15.056 

0.54 00:52:35.82 +56:37:33.6 17.766 1.674 14.829 13.875 13.289 



Table 5. The 2MASS, optical and X-ray counterparts of the IRAC sources searched within a match radius of 1 arc- 
sec. The radial distance is from the cluster center. The complete table is available in the electronic form only. 



Radial Distance 
(') 



a (2000) 
.m .s^ 



() (2000) 
to.ljl) 



3.6/xm 

(mag) 



4.5/im 

(mag) 



J 

(mag) 



H 

(mag) 



(mag) 



V 
(mag) 



(V-I) 
(mag) 



X-ray 
(Y/N) 



0.95 
1.00 
1.24 
1.24 
1.25 



00:52:33.90 
00:52:34.16 
00:52:31.18 
00:52:34.00 
00:52:31.32 



-56:37:11.8 
-56:37:04.2 
-56:37:16.5 
-56:36:46.3 
-56:37:12.1 



13.821 
12.265 
15.283 
14.390 
14.461 



13.391 
12.171 
15.177 
14.426 
14.463 



15.937 
13.391 
16.090 
15.233 
15.208 



14.981 
12.628 
15.485 
14.619 
14.566 



14.606 
12.366 
15.266 
14.460 
14.494 



19.876 
16.323 

17.617 
17.169 



2.169 
1.636 

1.409 
1.217 



N 
N 
N 
N 
N 



Table 6. Detail of identified cold IRAS point source. 



IRAS PSC 


RA (2000) 


Dec. (2000) 


-Pi 2 


F25 


^60 






^ h . rn . s ^ 


(o./.ll) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


00512+5617 


00:54:14.74 


+56:33:22.7 


1.98 


10.84 


44.81 


218.20 



o o o o o ; 



5" 



3 
c 



Tabic 8: MIR, NIR and optical data for Class 0/1, Class II and Class III objects as detected on the basis of MIR TCD (cf. 
§5). 



ID 


°2000 


^2000 


3.6 ^im 


4.8 it m 


5.8 ^im 


8.0 i±m 


J 


H 




V 


U — B 


B - V 


V - I 


X-ray 


IDfTablc 7) 




fh.rn.s-) 


fO.,.,1] 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Y/N 




Class 0/1 
































la 


00:52:10.99 


+ 56:30:58.8 


9.873 


7.856 


6.002 


4.992 


17.235 


15.499 


14.854 












213 


11. 


00:52:11.76 


+ 56:33:04.5 


10.409 


9.175 


8.023 


7.113 


17.067 


15.439 


13.445 












181 


Ic 


00:52:29.97 


+ 56:33:29.2 


11.129 


9.372 


8.114 


7.114 


















123 


Id 


00:52:17.15 


+ 56:33:42.5 


11.678 


9.538 


7.937 


6.846 


















137 


Ie 


00:52:26.58 


+ 56:33:25.5 


11.845 


10.677 


9.793 


9.170 


17.614 


17.246 


14.904 












127 


If 


00:52:24.56 


+ 56:33:50.1 


12.175 


10.345 


9.668 


9.714 


















121 


Ig 


00:52:13.48 


+ 56:33:41.4 


12.232 


11.233 


10.323 


9.459 


18.614 


16.133 


14.317 












156 


Hi 


00:52:16.45 


+ 56:31:45.9 


12.765 


12.256 


11.478 


10.362 


16.519 


14.960 


13.984 












200 


Ii 


00:52:15.92 


+ 56:33:49.5 


13.065 


11.618 


10.383 


9.644 


18.744 


16.025 


15.401 












139 


I.i 


00:52:19.56 


+ 56:32:58.4 


13.255 


12.723 


12.101 


10.955 


15.962 


14.941 


14.371 












16 1 


Ik 


00:52:32.54 


+ 56:32:33.6 


13.401 


12.822 


11.557 


9.870 


16.697 


15.463 


15.115 


19.743 






1.758 




1 19 


11 


00:52:21.37 


+ 56:30:40.4 


13.711 


12.544 


11.454 


10.546 


















209 


1m 


00:52:25.55 


+ 56:33:32.6 


15.221 


13.185 


11.631 


10.246 


















126 


In 


00:52:31.37 


+ 56:33:31.2 


15.625 


13.613 


12.470 


11.700 


















120 



Ha 


00 


52 


20 


84 


+ 56 


33 


OS 


2 


10 


048 


9 


305 


8 


535 


7 


695 


16 


288 


14 


617 


12 


502 


lib 


00 


52 


22 


39 


+56 


34 


26 


5 


10 


957 


10 


122 


9 


871 


8 


993 


14 


7 11 


13 


121 


12 


207 


lie 


00 


52 


12 


49 


+ 56 


34 


13 


7 


11 


148 


10 


449 


9 


714 


8 


798 


17 


305 


14 


889 


13 


359 


lid 


00 


52 


27 


34 


+56 


34 


03 


2 


11 


345 


10 


925 


10 


410 


9 


629 


15 


891 


14 


1 18 


12 


973 


He 


00 


52 


27 


31 


+ 56 


32 


04 


8 


11 


375 


10 


583 


9 


867 


9 


125 


15 


882 


14 


231 


13 


888 


Ilf 


00 


52 


35 


75 


+56 


34 


31 


1 


11 


391 


10 


820 


10 


330 


9 


896 


15 


383 


13 


871 


12 


827 


tig 


00 


52 


18 


83 


+ 56 


33 


05 


9 


11 


590 


11 


187 


10 


813 


10 


385 


14 


510 


13 


890 


12 


821 


Ilh 


00 


52 


19 


31 


+56 


31 


14 


1 


11 


814 


11 


269 


10 


714 


9 


96 7 


14 


582 


13 


179 


12 


947 


Hi 


00 


52 


21 


59 


+ 56 


31 


31 


8 


11 


947 


11 


517 


11 


108 


10 


127 


15 


1187 


13 


769 


13 


03 1 


II. i 


00 


52 


10 


89 


+56 


34 


06 


3 


12 


170 


11 


634 


11 


212 


10 


561 


17 


828 


15 


770 


1 1 


193 


Ilk 


00 


52 


33 


69 


+ 56 


35 


59 


8 


13 


010 


12 


775 


12 


496 


11 


841 


14 


976 


14 


000 


13 


657 


III 


00 


5 2 


19 


49 


+ 56 


32 


52 


7 


13 


169 


12 


694 


12 


028 


11 


484 


16 


189 


15 


358 


14 


957 


1 1 in 


00 


52 


1 1 


89 


+56 


31 


44 


8 


13 


885 


13 


340 


12 


702 


11 


638 


17 


1 17 


15 


553 


15 


261 



1 IV 
101 

1 11 

102 
170 

71 
159 
20 1 
198 
153 

29 
170 
202 



00 00 00 H O E 

tot -j k £t w » 

if. Ut W 00 Ul -■ 



*. to 03 O (0 I 
h oo u ui *i (r 

00 10 W 03 CO «- 



CO 

CL O 

O CD 



•2, 

3 

Q 

a 
D 



O 

o 

to 



o' 



ass III 
Ilia 


00 


52 


16 


82 


+ 56 


31 


33 


4 


8 


601 


8 


477 


8 


431 


8 


486 


8 


736 


8 


556 


8 


515 












Y 




mi- 


00 


52 


10 


31 


+56 


31 


31 


6 


9 


252 


8 


835 


8 


513 


8 


192 


10 


904 


10 


402 


9 


974 


12.828 


0.027 





893 


1.362 




207 


me 


00 


52 


14 


65 


+ 56 


34 


39 


6 


9 


699 


9 


79 1 


9 


63 1 


9 


569 


10 


558 


10 


000 


9 


839 


12.682 


0.790 


1 


145 


1.268 






Hid 


00 


52 


3 1 


20 


+56 


32 


21 





10 


173 


10 


23(1 


10 


200 


10 


208 


11 


137 


10 


1 19 


10 


287 


13.594 


1.090 


1 


352 


1.448 






I He 


00 


52 


36 


28 


+ 56 


36 


00 





10 


271 


10 


315 


10 


207 


9 


944 


10 


395 


10 


313 


10 


334 


10.867 


-0.507 





225 


0.362 






IHf 


00 


52 


42 


88 


+ 56 


35 


16 


3 


10 


712 


10 


727 


10 


642 


10 


789 


10 


997 


10 


835 


10 


768 


11.885 


0.240 





467 


0.572 






Illg 


00 


52 


10 


55 


+ 56 


35 


50 


1 


11 


22 1 


11 


230 


11 


111 


11 


205 


11 


859 


11 


548 


11 


391 


13.517 


0.549 





658 


0.933 


Y 




Illh 


00 


52 


25 


71 


+ 56 


3 1 


15 


5 


11 


913 


11 


300 


10 


783 


10 


523 


17 


138 


14 


700 


13 


611 












Y 


95 



o 

Q 

a 



O 
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Table 9. The mass and age of the YSOs having optical counterparts along with the associated errors. The ID is as same as in Table 7. 



ID 


Mass ±cr 


Age ±er 


ID 




Mass ±er 


Age ±cr 


ID 


Mass ±<7 


Age ±er 




(M 


0) 


(Myr) 






(M Q ) 


(Myr) 




(M 


0) 


(Myr) 


1 


2.20 ±0 


11 


1.01±0.21 


48 


1 


15 ±0.07 


0.58 ±0 


09 


192 


0.35 ±0 


02 


0.37 ±0 


27 


4 


1.11 ± 


08 


1.36 ±0.30 


50 





52 ±0.03 


0.72 ±0 


06 


197 


0.53 ±0 


04 


0.42 ±0 


1G 


6 


0.84 ±0 


05 


0.68 ±0.09 


51 





57 ±0.05 


1.89 ±0 


51 


199 


0.96 ±0 


07 


0.34 ±0 


04 


7 


2.72 ±0 


13 


3.09 ±0.39 


53 





89 ±0.07 


2.19 ±0 


55 


201 


0.45 ±0 


03 


1.01 ±0 


12 


8 


0.53 ±0 


05 


2.10±0.61 


59 





36 ±0.01 


0.31 ±0 


23 


206 


0.34 ±0 


01 


0.10±0 


01 


9 


6.29 ±0 


30 


0.36 ±0.04 


60 


1 


08 ±0.06 


0.43 ±0 


05 


209 


5.34 ±0 


07 


0.17±0 


02 


12 


1.84 ±0 


10 


1.19 ±0.24 


63 





54 ±0.04 


1.45 ±0 


34 


214 


0.81±0 


06 


3.10±0 


71 


14 


0.72 ±0 


07 


0.67±0.10 


65 





82 ±0.04 


0.92 ±0 


12 


217 


0.44 ±0 


01 


0.10±0 


00 


18 


0.97 ±0 


08 


3.71 ±0.97 


67 





91 ±0.06 


3.96 ±0 


82 


220 


1.48 ±0 


03 


4.86 ±0 


29 


22 


0.65 ±0 


05 


2.20 ±0.58 


70 


1 


87 ±0.05 


2.79 ±0 


57 


221 


1.70 ±0 


04 


3.82 ±0 


G8 


23 


0.73 ±0 


07 


1.69 ±0.40 


75 





53 ±0.03 


1.01 ±0 


14 


227 


1.43 ±0 


03 


4.85 ±0 


27 


25 


3.31 ±0 


14 


1.66 ±0.27 


76 


1 


09 ±0.06 


0.77 ±0 


12 


228 


0.94 ±0 


03 


5.00 ±0 


00 


27 


1.66 ±0 


08 


1.91 ±0.40 


79 





79 ±0.03 


5.00 ±0 


00 


229 


3.10±0 


09 


0.96 ±0 


25 


29 


0.75 ±0 


05 


0.77 ±0.07 


80 


1 


11±0.07 


2.67 ±0 


62 


230 


0.59 ±0 


04 


1.38 ±0 


27 


30 


0.84 ±0 


06 


1.59 ±0.34 


85 


1 


14 ±0.03 


4.99 ±0 


03 


231 


0.82 ±0 


04 


0.41 ±0 


03 


31 


0.49 ± 


03 


0.87 ±0.09 


87 


1 


09 ±0.07 


2.37 ±0 


53 


233 


1.22±0 


07 


0.39 ±0 


04 


32 


0.92 ±0 


06 


0.85 ±0.15 


88 





73 ±0.05 


1.24±0 


18 


234 


0.87 ±0 


03 


5.00 ±0 


00 


34 


0.51 ±0 


03 


0.78 ±0.07 


90 





87 ±0.05 


1.67±0 


34 


235 


0.64 ±0 


04 


4.48 ±0 


72 


36 


0.60 ±0 


04 


0.80 ±0.15 


96 


1 


39 ±0.05 


2.74 ±0 


55 


236 


0.52 ±0 


02 


0.10±0 


00 


38 


1.64 ±0 


04 


3.41 ±0.65 


108 





79 ±0.05 


1.68 ±0 


29 


242 


0.40 ±0 


02 


0.91 ±0 


08 


40 


0.49 ± 


01 


0.10±0.01 


125 





88 ±0.05 


0.85 ±0 


14 


245 


0.49 ±0 


04 


1.12 ± 


18 


43 


0.81 ±0 


04 


0.95 ±0.13 


149 





78 ±0.03 


5.00 ±0 


01 


250 


0.35 ±0 


01 


0.30 ±0 


23 


44 


1.34 ±0 


03 


4.96 ±0.05 


150 





78 ±0.03 


5.00 ±0 


01 


259 


1.03±0 


06 


0.97 ±0 


1G 


45 


0.90 ±0 


05 


4.39 ±0.72 


169 


1 


04 ±0.03 


5.00 ±0 


00 


268 


0.82 ±0 


03 


5.00 ±0 


00 


4G 


0.73 ±0 


05 


1.61 ±0.27 


175 





30 ±0.01 


0.12±0 


06 












47 


0.48 ± 


03 


0.57±0.11 


187 


1 


32 ±0.03 


4.97 ±0 


04 













Table 10. Statistics of probable CTTSs in three sub-regions. Numbers given in parentheses are in percentage. 



Region 


Total 


NIR-cxcess 


Ha 


Probable CTTSs 


A(H-K) 


A v 




stars 


stars 


stars 


(NIR-excess 


(mag) 


(mag) 










+ Ha stars) 






stars with error less than 0.1 mag 


1 


168 


2 (2.5) 


3 (3.8) 


5 (6.3) 


0.06 ±0.03 


1.6±0.2 


2 


107 


6 (33.3) 


1 (5.5) 


7 (38.9) 


0.08 ±0.04 


3.9 ±1.9 


3 


97 


4 (50.0) 


3 (37.5) 


6 (75.0)° 


0.16 ±0.05 


3.6 ±2.3 


Field 


89 












all stars independent of errors 


1 


425 


8 (3.3) 


3 (1.3) 


11 (4.6) 


0.12±0.08 


2.5 ±1.8 


2 


313 


32 (25.2) 


1 (0.7) 


33 (26.0) 


0.10±0.07 


4.4 ±2.8 


3 


241 


12 (21.8) 


3 (5.5) 


14 (25.5) a 


0.11±0.07 


4.6 ±2.8 


Field 


186 













a: One Ha stars have NIR-excess 



Table 11. CTTS fraction fcTTS = N C tts/( n CTTS +N wtts ). 



Region CTTS 




WTTS 




fcTTS 




error < 0.1 


all 


error < 0.1 


all 


error < 0.1 


all 


1 2 


8 


38 


56 


0.05 


0.12 


2 6 


32 


18 


39 


0.25 


0.45 
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Fig. 1. V band image of the NGC 281 region taken with the Kiso Schmidt. The big and small 

white boxes are the areas covered by the Ho and Spitzer observations and the black box repre- 
sents the Chandra observations respectively. The X and Y axes are in RA and Dec. in J2000. 




Fig. 2. (left panel) Isodensity contours of the 
tours are plotted above 3 sigma levels with the 
resents 17 stars/pc 2 . (right panel) Isodensity 



distribution of the 2MASS sources. The con- 

step size of 5 stars/pc 2 . The lowest contour rep- 

contours of the distribution of the identified YSOs. 




Fig. 3. Radial density profile of the cluster using the 2MASS data. The continuous curve show the least-square fit of the King 
(1962) profile to the observed data points. The error bars represent ±\/jV errors. The dashed line indicate the density of field stars. 
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Fig. 5. (V — K) vs. (B — V) TCDs for the nearby reference region (left panel) and for the cluster region (r < R c i) (right panel. Open 
and filled circles represent probable cluster members and field stars with normal reddening, respectively. The continuous lines show the 
least square fits to the distributions of the probable cluster members. The dashed lines show the fits to the distribution of field stars. 
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Fig. 6. Left panel: NIR TCD of the X-ray sources (open circles), Ha emission (star symbols), NIR-excess sources (open trian- 
gles), probable WTTSs (filled circles) and O-type star (open square), having photometric errors less than 0.1 mag, in the NGC 
281 region. Right panel: same as left panel but for all the sources in the reference region detected in the JHK S bands with 
the photometric errors less than 0.1 mag. The sequences for dwarfs (solid curve) and giants (thick dashed curve) are taken 
from Bessell & Brett (1988). The dotted line represents the loci of unreddened T Tauri stars (Meyer et al. 1997). Dashed 
straight lines represent the reddening vectors (see the text). The crosses on the dashed lines are separated by Ay = 5 mag. 
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Fig. 7. (left) IRAC MIR TCD of the detected sources. The sources lying within the box are Class II 
sources. The sources located around [5.8]-[8.0] ~ and [3.6]-[4.5] ~ are field/Class III stars. The sources 
with [3.6]-[4.5] > 0.8 and/or [5.8]-[8.0] > 1.1 represent Class 0/1 objects. The horizontal continuous line shows 
the adopted division between Class I and Class I /II sources (see Megeath ct al. 2004). (right) IRAC CMD 
for sources detected in the 3.6 and 4.5 fim bands only. Encircled sources represent objects with X-ray emission. 
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Fig. 8. The sample SEDs obtained by using the optical, NIR and MIR observations. Stars Ik, Ilk and Illg represent Class 
0/1, Class II and Class III sources, respectively, classified according to the MIR TCDs (cf. §4.3.2 ) as given in Table 7. 




Fig. 9. J I (J — H) CMD of the X-ray sources (open circles), Ha emission stars (star symbols), NIR-excess sources (open triangles), 
probable WTTSs (filled circles) and O-type star (open square), having photometric errors less than 0.1 mag, in the NGC 281 region. 
The isochrone of 2 Myr (Z=0.02) and PMS isochrones of age 1, 5 and 10 Myr by Marigo et al. (2008) and Siess et al. (2000), 
respectively, corrected for a distance of 2.81 kpc and reddening E(B — V) m i n = 0.32 mag are also shown. The parallel slanting 
dashed lines denote loci of 1 Myr old PMS stars having masses in the range of 0.1 to 3.5 Mq taken from Siess et al. (2000). 
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Fig. 10. V/(V— I) CMD for (a) stars in the cluster region and (b) stars in the reference region, (c) is a statistically cleaned CMD. 
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Fig. 11. Statistically cleaned Vo/(V — I)o CMD for stars lying in the cluster region (left panel). Filled circles (Ages < 5 
Myr) are used to estimate the MF of the region. Right panel shows the Vo/(V — I)o CMD for the Ha emission stars (star 
symbols), NIR-excess stars (triangles) and X-ray sources in the 'F' region of the NIR TCD (solid circles) (see §4.3 for de- 
tails) collectively. The ZAMS by Marigo et al. (2008) and the PMS isochrones of 1.2,5,10 Myr along with evolutionary 
tracks for different mass by Siess et al. (2000) are also shown. All the curves are corrected for a distance of 2.81 kpc. 
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Fig. 12. Histogram showing the distribution of ages of the YSOs. The inset shows the distribution of 

random errors. The error bar on Y-axis represent amount of scatter in each bin estimated on the ba- 
sis of errors associated with age estimates as given in Table 9 and using Monte Carlo simulations. 
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Fig. 13. A plot of the mass function in the cluster. log <f> represents log (JV/dlog m). The error bars represent the 
errors. The solid line shows the least square fit to the mass range 2.0 < M/Mq < 54 (shown by filled circles). 




Fig. 14. (a) Comparison of the observed KLF in the reference field and the simulated KLF from the star count modeling. 
The filled circles denote the observed iC-band star counts in the reference region, and the triangles represent the simula- 
tion from the Galactic model (see the text). The error bars represents the ±\/]V errors. The KLF slope (a, see §4.5) of 
the reference region (solid line) is 0.36 ± 0.05. The simulated model is also gives a similar value of slope (0.34 ±0.02). (b) 
The corrected KLF for the probable members in the cluster (see the text). The straight line is the least-square fit to the 
data points in the magnitude range 10.5-14.25. The KLF slope (a) for the cluster region (R c i < 5 arcmin) is 0.37 ±0.07. 
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Fig. 15. The IRAS intensity maps for the cluster region in 12 fim (top left), 25 /im (top right), 60 fim (bottom left) 
and 100 fim (bottom right). The contours arc at 20, 25, 30, 35, 40, 50, 60, 70, 80 and 90 % of the peak value 
of 20 MJy/ster, 48 MJy/ster, 412 MJy/ster and 907 MJy/ster in 12, 25, 60 and 100 fim respectively. The loca- 
tions of IR-excess stars (probable CTTSs, triangles), Ha emission stars (star symbols), IRAS point sources (crosses), 
and the O-type star (open square) are also shown in the image. The abscissa and the ordinates are in the J2000. 
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Fig. 16. The locations of IR-excess stars (probable CTTSs, triangles), Ha emission stars (star symbols), IRAS point sources 
(crosses), the O-type star (open square), and WTTS (filled circles) are overlaid on the DSS-2 Ft, band 30 X 30 arcmin 2 image. The clus- 
ter region is represented by the circle having the center "C". The CO contours taken from Hcnning et al. (1994) are shown by the blue 
dashed lines and the NVSS (1.4 GHz) radio continuum contours by the white lines. The radio contours are 5, 10, 15, 20, 30, 40, 60, 
80% of the peak value of 0.04 Jy/Beam. Three subrcgions are also marked by the boxes. The abscissa and the ordinate are in J2000. 
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Fig. 17. Variation of the age and NIR-excess [A(H — K)\ of the YSOs as a function of radial distance from the O 
star HD5005. The symbols are same as in Fig. 6. The dotted line represents the Mean ± a of the distribution. 
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Fig. 19. NIR TCDs of the YSO candidates in the three subregions. The sequences for dwarfs (solid curve) and giants (thick 
dashed curve) are taken from Bessell & Brett (1988). The dotted line represents the loci of unreddened T Tauri stars (Meyer et al. 
1997). Dashed straight lines represent the reddening vectors (see the text). The crosses on the dashed lines arc separated by Ay 
= 5 mag. Open squares show the data having optical counterparts used in Fig. 18. The other symbols are as same as in Fig. 6. 
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Fig. 20. Spatial distribution of the Class 0/1 (filled triangles) and Class 11 (open triangles) sources detected by Spitzer 
IRAC overlaid on the 2MASS image of the NGC 281 region. The square box represents the area observed by Spitzer 
IRAC. The dashed curves indicate the boundaries of the three sub-clusters named as a, b and c (see §5 for details). 
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Fig. 21. Spitzer [3.6]-[4.5] MIR colour vs. Chandra X-ray source median energy. Filled cir- 

cles and squares represent Class II and Class III sources, respectively.. The dashed lines indicate 

the boundaries of the distribution of similar objects in the IC 1396N region (Getman et al. 2007). 



